Coexisting Kondo singlet state with antiferromagnetic long-range order: A possible 

ground state for Kondo insulators 
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The ground-state phase diagram of a half-filled anisotropic Kondo lattice model is calculated within 
a mean-field theory. For small transverse exchange coupling J± < Jxci, the ground state shows an 
antiferromagnetic long-range order with finite staggered magnetizations of both localized spins and 
conduction electrons. When J± > J±c2, the long-range order is destroyed and the system is in 
a disordered Kondo singlet state with a hybridization gap. Both ground states can describe the 
low-temperature phases of Kondo insulating compounds. Between these two distinct phases, there 
may be a coexistent regime as a result of the balance between local Kondo screening and magnetic 
interactions. 

PACS numberes: 71.28.+d, 72.15.Qm, 75.20.Hr 



The Kondo lattice model is often considered as a the- 
oretical model for heavy fermion materials. For this 
model, an important issue arises from the interplay be- 
tween the Kondo screening and the magnetic interactions 
among localized spins mediated by the conduction elec- 
trons. The former effect favors a nonmagnetic Kondo 
singlet state in strong coupling limit, while the latter 
interactions tend to stabilize a magnetically long-range 
ordered state in weak coupling limit. The nature of such 
a transition between these two distinct phases has been a 
long standing issue since it was first suggested by Doniach 
0. The one-dimensional model was intensively studied 
at half filling, showing that its ground state is always a 
disordered Kondo singlet state (or a spin liquid state) 
1^. In higher dimensions, however, both antiferromag- 
netic long-range order (AFLRO) and disordered Kondo 
singlet states may occur Q. Recently, there have been 
more indications of such a transition from various ap- 
proximate treatments for the Kondo or Anderson lattice 
models, including variational Monte Carlo calculation [Q, 
higher-order series expansions ||] , quantum Monte Carlo 
simulations &|7i, and infinite dimensional calculations 



In this paper, we would like to consider the issue 
whether the disordered Kondo singlet state can coex- 
ist with an AFLRO at half-filling. First of all, we in- 
troduce an anisotropic Kondo lattice model by distin- 
guishing the longitudinal spin exchange interaction from 
the transverse one, because we notice that the longitudi- 
nal interaction describes a polarization of the conduction 
electrons by the localized spins, while the transverse one 
describes a spin-flip scattering of the conduction elec- 
trons off the localized spins. The former interaction is 
the origin of the magnetic interactions among localized 
spins, leading to an AFLRO; the latter interaction is re- 
sponsible for the local Kondo screening effects, yielding a 
disordered Kondo singlet state. We also notice that the 
AFLRO and local Kondo singlet order operators form an 
irreducible representation of an SU(2) algebra, which can 



be regarded as the spectrum generating algebra Q of the 
half-filled Kondo lattice model. Within the framework of 
a mean field theory, the magnetic interactions and the 
Kondo screening are considered on an equal footing, and 
the ground-state phase diagram of the model Hamilto- 
nian is calculated. 

The symmetric Kondo lattice model with anisotropic 
exchange couplings is defined as: 

H = J2 «^kcLck^ + ^11 + H± 

kCT 

Jl 

1 ^ 



(1) 



where the pseudo-fcrmion representation for the localized 
spins 5f = {dl^d,^ -dl^d,i)/2, S~ = dl^d,^, 5+ = dl^d^, 

has been used with a local constraint dj^dii +dl^dii = 1. 
Note that H± can also be rewritten in the form of 



-Jl 



describes the polarization of conduction electrons by 
the local impurity spins, leading to a magnetic instability, 
while H± corresponds to the spin-fiip scatterings, giving 
rise to the local Kondo screening effect. The latter effect 
has been investigated by various approaches, in particu- 
lar, those based on a 1/N expansion (N is the degen- 
eracy of the localized spin). However, the former effect is 
not treated on an equal footing in these approaches be- 
cause the magnetic interaction occurs there only at 1 /N^ 
order [ pT[ |. 

It has been known that the spectrum generating al- 
gebra plays an important role in analyzing complete 
spectra for a model. For instance, it has been used to 
study collective excitations and phase transitions in one- 
dimensional metals and two-dimensional Hubbard model 
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Here the spectrum generating algebra for the half- 
filled Kondo lattice model is given by generators of an 
SU(2) Lie group: 



i 

r"^ = i ^ [(4tC»T - cliCji) - {d\^d^^ - dj^dii) , (2) 

i 

which satisfy the commutation relations: [r^ , t^] = ±t^ , 
[t+, r~] = 2t^, the SU(2) algebra. Q is the AF reciprocal 
vector. We also find that an irreducible representation 
of this SU(2) algebra can serve as the order parameter 
operators for this model. They are given by: 

i 

i 
i 

with the following commutation relations [t^,K^] = 
2M^ [r±,M^] = K"^, [r^ii:±] = iif^, and [r^ M^] = 
0, where are the order parameter operators describ- 
ing the local Kondo singlet state and is the staggered 
magnetization operator describing the AFLRO in terms 
of a commensurate spin-density wave. The expectation 
values of these order operators are the corresponding or- 
der parameters. Note that the model Hamiltonian does 
not satisfy this SU(2) symmetry defined by the spectrum 
generating algebra, which can however provide us with a 
guide to choose the order parameters in the model Hamil- 
tonian. 

Now we introduce a mean-field theory treating equally 
the magnetic interactions between the localized spins and 
the itinerant electron screening aspects. The longitudinal 
interaction term is approximated as: 
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where the staggered magnetizations have been intro- 
duced as the AF order parameters mc ~ ~\ < — 

c\^c^l > e'^^ "' and rud ^ \ < d\^di^ - > e-'*^'""', 

and a minus sign has been absorbed into the definition 
of ruc for convenience. The transverse interaction term 
is approximated to be: 



Ji V 



with the hybridization order parameter describing the 
local Kondo singlet state —V =< cl-^di^ + dl^Cii > 
—< d-|Ci| + cl^dii > . Moreover, we could also intro- 
duce a d-electron chemical potential term Ej.o- Edd\^dic 
to fix the d-electron density to be one on each site, but 
actually this is not necessary at half filling because the 
electron-hole symmetry automatically imposes Ed = 
[^. Therefore, the mean field Hamiltonian can be writ- 
ten in the following form 



\mcmd 



Ji V 



), (3) 



where ^'[^ = (cL, 4+q<t)' transposition 

^'ko-, arid the matrix 



-fmda, -Ck, 
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(4) 



The static staggered magnetizations partially break the 
translational symmetry and the Brillouin zone is folded in 
half, so that the summation over k is taken in the reduced 
Brillouin zone. The quasiparticle bands are determined 
by the equation: | EI — H \— giving rise to four bands 
with dispersions. 



E±±{\^) = ± 
with 

i?'(k) = 



1 



J2 {ml + m2 )/4 + J2 y2 /2 ± e'{]s) 



ti + J({mi + m'^)/A + JiV'/2 

2 



t2 2 2 



and the ground-state energy is given by 



1/2 



Eg^2Y,' [E—(y) + i?-+(k)] +Af{J\\m,md + ^^), 

k 

corresponding to completely filling the two negative en- 
ergy bands with electrons. The self-consistent equations 
are obtained by minimizing the ground-state energy Eg 
with respect to m^, rric, and V , respectively. After the 
summations over momenta are transformed into integrals 
over energies by assuming a constant density of states of 
the conduction electrons in [— the self-consistent 
equations are expressed in the form 

A|| \ 2F2{e)md + {X\mcmd + X\V'^)mc ^ 



2Fi{e)F2{e) 



All f 2{F2{e) + 2e^)mc + {Xjmcmd + \'iV'^)md 

de s = rud, 



2Fi{e)F2{e) 
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r 2F2{e) + {Xlm,ma + XlV^) _ 
-1 



(5) 



where we have used the simphfied expressions in terms 
of Fi and F2, 

Fi{e) = ^e^ + A|(to2 + ^2^/4 ^ xlvy2 + F^ie), 
F2{e) = (\\mcmd + \\V^)\xlmle\ (6) 

and we have also used the notations: A[| = J\\ID, 
Aj^ = J±/D, where 2D is the bandwidth of the con- 
duction electrons. 

When A|| ^ Ai^, the polarization effect described by 
H\\ dominates, and the ground state has an AFLRO char- 
acterized by finite staggered magnetizations for both lo- 
calized spins and the conduction electrons. Thus, the to- 
tal staggered magnetization of the system should be given 
by jM = (TO(J -I- ?7ic), while the staggered magnetic 
moment is /i = (m^ — m^. Let us assume for the mo- 
ment they are the only order parameters, i.e. Y = 0. All 
equations involved are greatly simplified. Due to unit cell 
doubling the conduction electron band is folded in half 

to form two new bands: e±(k) = 4- {.J\\md)^ jA,. 

There is a gap in the conduction electron excitation spec- 
trum: A = J\f^md/2. The charge gap is twice this value, 
Ach = J\\md- Moreover, the ground-state energy of the 
AFLRO state can also be found 
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(7) 



where M is the total number of the lattice sites, and 

In the opposite limit, when A^ S> A||, the Kondo 
screening effect described by dominates, and the 
ground state has a finite local Kondo order parameter 
V . If we assume again for the moment it is the only or- 
der parameter, the system is a nonmagnetic band insu- 
lator with an effective hybridization: = A^ sinh (^xt) • 
A similar result has already been given by the 1/N ex- 
pansion approach [|l^, which is believed to describe the 
correct low-energy Kondo physics. The quasiparticle 
excitation spectrum is expressed in the form e'-t(k) = 

i ^Ek ± \/^k + {J^yT^ 1 ^i^d there is a small hybridiza- 
tion gap Ahy — {J±V)'^ /2D, which splits the Kondo res- 
onance formed at the Fermi level. At half filling, each 
lattice site has one conduction electron and one local- 
ized spin, and they can form a local singlet state, thus 
the system becomes a disordered nonmagnetic insulator, 
a band- insulator with both charge and spin gap 2Ahy 
The ground-state energy can be also calculated 



-f- = --coth( — 
AfD 2 V A 



(8) 



Physically, this disordered Kondo singlet state is adia- 
batically connected to the usual Kondo spin liquid state 

0- . . 

Comparing the respective ground-state energies 
and , we find the phase boundary between the two 
ground states: 



coth 




, (9) 



which is displayed by the solid line in Fig.l. Below this 
line, the AFLRO state is more stable E^^ < E^^; above 
this line the disordered Kondo singlet state is more stable. 
We note that this line intersects the diagonal Ax = A|| 
at 0.58. Taking into account that the cut-off parameter 
D = 2t with t as the nearest neighbor hopping, in two 
dimensions, where our assumption of a constant density 
of states is better justified, we find Jc/t — 1.16. This 
value is not far from the numerical results Jc/t = 

1.40 ^ 1.45, obtained for the isotropic model. 

Now we turn to discuss the main new result of this 
paper, namely the possible coexistence of AFLRO with 
Kondo singlet state. In fact, our system of self-consistent 
equations allows solutions with all order parameters 
TOd,TOc and V being nonzero. Let's start from the re- 
gion where AF order dominates, with finite rud and rric 
to see where the instability towards the Kondo singlet 
state emerges. Assume V being small, but non-zero, we 
obtain the following solution. 




^tanhf^A 

1 + ^ \md\\\ 



, (10) 



where 



This instability 



rrid ~ 4 

line for the AF state is delineated by the thick dotted line 
in Fig.l. In a similar way, we can also determine the in- 
stability boundary in the Kondo singlet phase, which cor- 
responds to the appearance of small AF order parameters 
TOc and TTid- From the above self-consistent equations, we 
can derive the critical value of the coupling parameter 



A I sinh I ^ 



4Ai 



A^sinhf^") -1-1 



(11) 



This 



and a relation between m„ and m^: — 

instability line for the disordered Kondo singlet state is 
displayed by a thin dotted Hne in Fig. 1. 
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Between these two instability lines (Aj_ ci < < 
A_L c2), there is a narrow regime where the AF and the lo- 
cal Kondo singlet screening order parameters coexist to 
balance the energy gain from the transverse and longi- 
tudinal exchange couplings. The disordered Kondo sin- 
glet state and the AFLRO state can both be used to 
describe two distinct ground states of the conventional 
Kondo insulators. It seems to us that a new phase may 
be present when the exchange coupling parameters are 
tuned carefully. In such a new phase, the dynamical mag- 
netic structure factor should have two contributions: a q- 
indcpendent single site slow component, which is typical 
of the localized Kondo-type excitations, and a strongly 
q-dependent intersite fast component, reflecting the mag- 
netic interactions. These features can be detected in in- 
elastic neutron scattering experiments for some Kondo 
insulating materials. Moreover, since the localized spins 
are partially screened by the conduction electrons, and 
the residual localized spins still have weak AF long-range 
correlations mediated by the polarization of the conduc- 
tion electrons. This novel feature may be related to the 
small magnitude of the induced staggered magnetic mo- 
ments for URu2Si2 and UPts, the so-called heavy fermion 
micromagnetism [ p^ . 

We realize that our result is based on a mean field cal- 
culation, and its validity should be further checked by 
more rigorous analytic and numerical treatments. Even- 
tually, it should be verified by experiments. We also 
note that the numerical results for the isotropic model 
indicate a second order phase transition between 
the Kondo singlet and AFLRO states. However, the or- 
der of phase transitions is a very sensitive issue, and the 
isotropic case may well be a degenerate point of the more 
general anisotropic model we consider here. 

In conclusion, we have considered a half-filled anisotropic 
Kondo lattice model within a mean field theory. The 
ground-state phase diagram has been calculated. In ad- 
dition to the AFLRO phase and the Kondo singlet phase, 
we have found that both of these two distinct phases can 
coexist as a result of the balance between the Kondo 
screening effects and the magnetic interactions, which 
provides a possible new ground state for the Kondo in- 
sulating compounds. 
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Figure Caption 

Fig.l. The ground-state phase diagram of the half- 
filled anisotropic Kondo lattice model in the — A|| 
plane. The narrow area between the lines of (A_l)ci and 
(A_l)c2 is the regime where the local Kondo singlet state 
and the AFLRO may coexist. The solid line corresponds 



to the boundary of i?, 



AF 



a ' 



and the thin dashed 



line denotes the isotropic limit of the model. 
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Figure 1 
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